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Abstract 



We present the results of the complete one-loop computation of the B^^ ^ l^l^ decay rate 
in the MSSM. Both sources of the FCNC, the CKM matrix and off-diagonal entries of the 
sfermion mass matrices are considered. Strong enhancement of the branching ratio (compared 
to the SM prediction) can be obtained in the large tan/? ~ rrit/mh regime in which the neutral 
Higgs boson "penguin" diagrams dominate. We make explicit the strong dependence of this 
enhancement on the top squarks mixing angle in the case of the chargino contribution and 
on the fi parameter in the case of the gluino contribution. We show that, in some regions of 
the MSSM parametre space, the branching ratio for this process can be as large as lO"*^^"^) 
respecting all existing constraints, including the CLEO measurement of BR{B — »• Xsj). We 
also estimate, that for chargino and stop masses ~ 0(100 GeV) BR{B^ — > l^l'^) with W = er 
or fiT can be of the order of 10~^^ for the still allowed values of the off-diagonal entries in the 
slepton mass matrix. 
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1. Introduction 



Extensions of tlie Standard Model (SM) usually predict new contributions to the flavour 
changing neutral current (FCNC) processes. For example, adding in the most general way a 
second doublet of the Higgs flelds to the standard theory of electroweak interactions typically 
leads to large amplitudes of FCNC processes mediated at the tree level by neutral Higgs par- 
ticles. Restricting appropriately the possible form of couplings of the two Higgs doublets to 
up- and down-type fermions eliminates such tree level contributions to FCNC processes but, of 
course, new contributions induced by loops involving the physical charged scalar still remain. 
Charged Higgs boson contributions to FCNC processes depend however on the same elements of 
the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix as does the standard boson con- 
tribution and, thus, amplify the effects of the FCNC source that is present in the SM, rather 
than being an independent new source of such processes. Nevertheless, requiring the effects 
of the charged Higgs boson not to spoil succesful predictions of the standard theory leads to 
interesting bounds on the {Mh+, ta.n (3) plane where tan/5 = V2/V1 is the ratio of the vacuum 
expectation values of the two Higgs doublets. In particular, in the popular two Higgs doublet 
model of type II (2HDMII), in which the first doublet couples only to leptons and down-type 
quarks and the second one couples to up-type quarks only, processes like B Xg'-f and K^K^ 
mixing together with bb constrain the plane (M/^+, tan/5). In particular, for tan (3 > 3 

B Xgj requires if^ to be heavier than ~ 165 — 200 GeV 

Supersymmetric models like the MSSM, which are the most popular and best motivated 
extensions of the SM, apart from containing a charged Higgs boson H^, induce yet additional 
contributions to the FCNC processes. Firstly, in such models the effects of the CKM mixing 
can be further amplified through the loops involving charginos and squarks. Secondly, as there 
is no reason why the squark mass matrices should be diagonal in the same (so-called super- 
CKM) basis as quarks, the sfermion sector of such models is in general a new, independent of 
the CKM matrix, source of the FCNC processes. 

Current experimental data on FCNC processes provide important constraints on these 
sources of fiavour nonconservation in supersymmetric models. (Extensive reviews are the refs. 
and Taking the CKM matrix as the only source of the FCNC processes, the current ex- 
primental data on B ^ Xg'y and K^K^, B^B^ mixings impose some constraints on the MSSM 
parameter space. These constraints, which correlate masses and composition of charginos and 
top squarks with the mass of the charged Higgs boson, depend in part on the element Vtd 
of the CKM matrix which is not directly measured 0, ^ and become weaker with growing 
sparticle masses. Effects of the nonzero off-diagonal entries of the sfermion mass matrices are 
usually analyzed separately 0, ^ . Stringent constraints apply to the entries causing transitions 
between the first two generations. Bounds on the entries connected to the third generation, 
which follow from the B^B^ mixing and B Xs'~f decay are much weaker. Thus, large devia- 
tions from the rates predicted in the SM can still be discovered in the forthcoming (or already 
running) experiments like BaBar (SLAC), BELLE (KEK), CLEO (Cornell), HERA-B (DESY) 
and LHC (CERN). In this context, particularity interesting process to look at are the decays 
^s(d) ~^ because they are clean theoretically being almost free of hadronic uncertainties. 
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Several papers analysed this process in the MSSM 0, H, 0| under various assumptions and 
with different approximations. In this paper we perform a complete calculation Q of the process 
— ^ in the MSSM with emphasise on qualitative understanding of the dominant effects. 
To this end we derive simple analytical formulae approximating the main contributions. We 
reconfirm that for values of tan/? < 20 the rate of this process is not significantly enhanced 
compared to the prediction of the SM (apart from the case of tan/3 ~ 0.5 and light if^ 
which is not favoured theoretically within the supersymmetric framework). In agreement with 
earlier papers [§] we find that large enhancement of the branching ratio is obtained in the case 
of large tan (3 values due the neutral Higgs boson penguin graphs. This has been previously 
made explicit in ref. (possible role of such contributions to K^-K^ and B^-B^ mixings has 
been emphasized in ref. [^) in which the contributions of charginos as a source of the flavour 
changing has been considered. We demonstrate strong dependence of the decay rate on the 
value of the stop mixing angle and explain it using our analytic formulae. Moreover, we extend 
previous calculations by analysing also the case of the flavour mixing induced by squark mass 
matrices. In the latter case we find very strong dependence on the /i parameter. Finally we 
correlate the predictions for B^^ ^ with the constraints imposed on the parameter space 
by other processes, in particular by the measurement by CLEO of the B — > Xg'y branching 
ratio. We find that even respecting all those constraints, BR{B^ — /i^/x^) can be enhanced 
up to 10~*^^~^) for tan/3 ~ mt/mt, ^ 50. Moreover, for such values of tan/?, and the off- 
diagonal 13 entries of the down-type squark mass matrix saturating the existing bound p, |[, 
also BR{B^ /^~/^^) can be of the same order of magnitude. This means that the unsuccesful 
search done at CLEO [jlT| already provides a constraint on the 13 off-diagonal entries of the 
down-type squarks which, for some values of the other MSSM parametres, is stronger than 
the one given in |^. Finally, we also estimate that for chargino and stop masses ~ (9(100 
GeV) B^ with //' = er or /ir can be of order 10~^^ for the still allowed values of the 

off-diagonal entries in the slepton mass matrix. 

2. General structure of the amplitude and the SM prediction 

The effective Lagrangian describing the didj — > Ia^b transition has the general form 

^eff = ^CxO^ (1) 

X 

in which Ox are the local four-fermion operators and Cx are their Wilson coefficients (we 
suppress quark and lepton flavour indices on Ox and Cx as well as on various formfactors which 
will appear in the following). Four vector- vector operators 0\y = {dj'j^Pxdi)(JB'y'^PYiA) 
and four scalar operators Oxy = {djPxdi){lBPYl a) (where X,Y = LL, RR, LR and RL) 
contribute to this process. In addition, two tensor operators exist but they do not contribute to 
this process (their matrix elements vanish when taken between one meson and vacuum states). 
In the following we will specify the formulae to the case of -8^(5) = bd{s) decay, hence we will 
take J = 3 and I = 1 = d for B'^ or I = 2 = s for B^. Furthemore, due to the pseudoscalar 

^In what follows we display only formulae for the dominant contributions. The plots arc, however, based on 
the programme including contributions from all relevant one- loop diagrams. 
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nature of the mesons we need only two matrix elements^: 

(0|67Vc^/|5?(g)) = -^fB,q^ 
{0\h'dj\B'j{q)) 



+ifBj 



Ml 



nidi + 



(2) 



Using one finds the total width 
^ Mb 



167r 



xl) {lap [l - {xA - xbY] + l&r [l - {XA + XB?] } 



(3) 



where f{x,y) = - 2{x + y) + {x - yf, xa 
in terms of the Wilson coefficients as 



itii^/Mb and the coefficients a and b are given 



Ibi ^ 
4 

Ibi 



{mi„ - mi. 



RR 



'RL 



Ml 
nib 



Ml 
nib 



'RR. 



RL 



/^b I /^b /^b 
^LL "T ^r.R ~ ^RR ~ ^ 



'LR 



RR 



'RL 



(4) 
(5) 



where we have neglected ni^j compared to nib- 



Three groups of diagrams contribute to the Wilson coefficients: Box diagrams, penguin 
diagrams and neutral Higgs boson penguin diagrams 0. Denoting the self energy diagrams on 
the external quark lines as —iT,{^), with 



(6) 



and vertex corrections to the couplings djdjZ^ , djdjS^ and djdjP^, where S'°(P°) is a neutral 
scalar (pseudoscalar) , respectively as 



+ ij>^{F^PL + F]^Pn) 
-z{F^PL + F^Pn) 
-{F[Pl + FEPr), 



(7) 



one finds (in the approximation S(p^) = S(0), F{q^) = F{0)) the following expressions for the 
Wilson coefficients generated by various penguin diagrams: 



C 



V 

XY 



V„e 



2swCwMz 



-F^r' 
2 ^ X^Y^ 



X,Y = L,R 



^The second follows from the first one by using the QCD equation of motion for the quark field operators; 
this fixes their relative sign. 

^At the one- loop level the photon penguin diagram does not contribute for the II final state due to the vector 
current conservation. As long as neutrinos are massless, the W final state can appear neither in the SM nor 
in the 2HDM; in the MSSM the II' final state can only be due to box contribution provided the slepton mass 
matrices remain non-diagonal in the lepton mass eigenstate basis. 
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from the penguin diagram, with c| = 1 — 2s^, c|j = — 2s^ and Sw (cvy) is the sine (cosine) 
of the Weinberg angle; 

Cf. = Cf« = £ cl« = Ci, = £ -i^l^LPI™, (9) 

from neutral scalar penguin diagrams; and 

1 '/'Ik 1 

cL = -cL = E M^T^^r-^ = -^1^ = E M^^^-'-^ (10) 

from neutral pseudoscalar penguin diagrams. We use here (and throughout) the notation of 
ref. Q in which = (/i°,if°), H^+^ = (A°,G°), = {H^,G^) and Z]f (Zf ) denotes the 
projection of the k-ih physical neutral CP-even (-odd) Higgs boson onto the real (imaginary) 
part of the neutral component of the Higgs doublet that couples to the down-type quarks. In 
addition, since at one loop penguin graphs cannot generate transitions —>■ W, we have set 
mi^ = mig = mi. In these formulae 

FlR = FlR + W7^ciR^lR (11) 



where 4 = 1- 2s^/3, 4 = -2s^/3, 



Fin = Fin - (12) 



and 



ylk ylk 

Ff = Ff + ^Sf, F^^ = F|-^E| (13) 

are the full effective vertices including the effects of flavour changing self energy diagrams on 
the external quark lines. Box diagram contributions to the Wilson coefficients can also be 
easily found. From eqs. (^,|5|,|8|-p!0D one sees that scalar penguin diagrams contribute only to h 
in eq (§) whereas the coefficient a receives contributions from both and the pseudoscalar 
penguin diagrams. The relative sign of the Z^ and the neutral Goldstone boson contributions 
to a should be such that the total contribution is independent of the gauge chosen for the Z° 
propagator. This is the case if 

- mjFl^ + miFl^ = -MzF[^n (14) 

for P referring to the Goldstone boson. Since the formfactor F^j^ for the physical pseudoscalar 
is related to the one for by the SUlC^) symmetry, this relation tests also the relative 
sign of the Z° and pseudoscalar penguin diagrams. 

The SM contribution to the II decay is well known [T^ (see also |]l4l)- The Higgs 



boson couplings to fermions are not enhanced so the scalar and pseudoscalar penguins are 
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negligibly small. The only important box diagram is the one with two W"^ which contributes 
only to CYl 



1 f -^f/ 



167r2 \swJ M2 4 



1 _^ log Xt 



Xt 



(15) 



where Xt = {nit/MwY and \ti = V^jVu. The effective djdjZ^ vertex receives contributions 
from loops involving both and the charged Goldstone bosons. One finds 



16n'^ 4:S^rCw 



Xt — 6 3xt + 2 



1- Xt (1 - Xtf 



\ogXt 



(16) 



and Fj^ = in the limit of m^^ = 0. Adding all one gets 



5i?(5° ^ 11) = T{Bf, 



2 r2 



TT 



\ |2 



nil 



(17) 



where t{B^) is the lifetime of the meson and |T^, |T4| 



Yoixt) 



Xt 



Xt 



Sxt 



Xt - 1 {xt- ly 



log Xt 



(18) 



In general, taking QCD corrections into account consists of computing corrections to the 
Wilson coefficients at the scale ~ Mz, and subsequently evolving the latter from the electroweak 
scale down to the hadronic scale Hh ~ n^'b- The first step of this procedure amounts to replacing 
Yo{xt) by Y{xt) = Yo{xt) + (as/47r)Fi(xt) ||TB|] where now Xt = {mt{mf)/MwY- Y{xt) can be 
conveniently parametrized as |l^ 



Y{xt) = 0.997 



166GeV 



1 1.55 



As far as the evolution is concerned, it has been noted in [|16l that the vector operators con- 
tributing to B^ II have zero anomalous dimensions. Hence, their Wilson coefficients do not 
evolve at all, whereas the evolution of the Wilson coefficients of the scalar operators result in 
multiplying them by mh{fih)/'m'b{Mz)- Consequently, if C^y are proportional to mb{Mz), their 
evolution is taken into account if this factor is replaced by mh{fih) which in turn cancels out 
with the factor l/mh{^ih) present in eqs. (^J^). As it will be apparent, whenever the coefficients 
CxY 3-^6 large, they are indeed proportional to mh{Mz). In the SM, including QCD corrections 
one finds 



BR{Bl /x/x) = 4.1 X 10" 



r{Bs) 
1.54 ps 



n 2 r 



245MeV 



Wt 



ts\ 



0.040 



mt{mt) 
166GeV 



n 3.12 



(19) 



3. Contribution of the extended Higgs sector 

As remarked in the introduction, the presence of the physical charged Higgs boson in the 
extended Higgs sector of the MSSM (or 2HDM) in general enhances the FCNC transition 
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rates generated by the CKM mixing matrix. This enhancement can appear through the 
contribution to box diagrams, penguin diagrams and neutral Higgs boson penguin diagrams. 
The latter type of diagrams can only be important in the large tan /3 > 30 regime in which the 
neutral Higgs boson couplings to the down-type quarks and charged leptons are enhanced by 
tan/3 factors. 

For low values of tan /9 < 20 the neutral Higgs boson penguin diagrams are small. It is also 
easy to check, that for such tan (3 values no box diagram can give significant contribution. Thus, 
the only large contribution can be due to the contribution to the penguin diagrams. 
Computing the relevant self energy diagrams (vector parts thereof) and vertex corrections one 
arrives at 



1 



167r2 {swCwY 



\ti cot^ (3 



1 yt 



Ml A I 



yt 



167r2 (swCw)^ 



Xti tan^ P 



mbiridj 1 yt 



Ml 



^i-yt 



1 + 



1 + 



1 



yt 



log yt 



l-yt 



logyt 



where yt = {rrit/Mfj+y. Taking into account only AF^ which is enhanced for tan/3 < 1 
amounts to replacing Y{xt) in eq. ([T7|) by 



Y{xt)~^Y{xt)-cot^P 



xt yt 



yt 



1 + 



i-yt 



logy* 



(20) 



The new contribution has the same sign as Y{xt) and, therefore, enhances the SM contribution. 
For example, for tan/3 = 0.5 and Mfj+ = Mw, BR{B^ II) is enhanced by a factor of 

1 + ~ 6.6 compared to the SM prediction. 



( 



For large tan /? ~ rrit/mb in the case of decay, F^, despite being suppressed by one power 
of nis/Mw, is about two orders of magnitude larger than and, for Mh+ ~ 100 GeV, is of 
the order of the SM contribution. However, in this regime, there are other contributions which 
are more important [|l^ 0. 



Firstly, the mixed, W^H^, box diagram in which couples to the 6-quark is also (9(tan^ /3) 



and is not suppressed by mdj/Mw 
gives |1^ 

1 / e \4 



After summation over different types of virtual quarks it 



C 



LR 



167r2 



SwJ 



mrrib 2^1 
A(/tan p- 



Xt 



Axh - xt 



xh -1 



log Xh - 



xt-l 



log Xt 



(21) 



where xh = {,Mh± / My/Y . The other W^H^ box is proportional to nidj/Mw and hence it is 
less important. The box diagrams containing two chaged scalars (either physical or Goldstone) 
are suppressed always by {mi/My/Y- Therefore, although the H^H^ box grows as tan^/?, it is 
not important even for tan (3 ~ 50. 



Secondly, there are neutral Higgs boson penguin diagrams. It turns out [0, that the 
dominant i.e. ~ tan/? part of the genuine djdiS^{P^) vertex correction cancels out and the 

^This cancellation is even simpler in the case of the MSSM than in the case of the 2HDM(II) considered in 
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only contribution arises from the scalar parts of the self energies of external quarks (for 
decay it is S;^ which is dominant): 



yS 



IGtt^ \sw 



rubXtiXt 



xh - Xt 



Xh 



Xh 



hgXH - 



Xt 



Xt 



log Xt 



(22) 



where xh = {Mjj± / M^) 



Using the formulae ([TT|.13) and the fact that in the MSSM for neutral CP-even scalars for 
large values of tan (3 the following relations hold^ 



sm a 



Ml ^ ...^ 

r2 i\ 4-2 



Ml for Ma < Mz 



cos^ a ^ 1, MfjK. Ma for Ma > Mz 



(23) 



= Ml + M, 



w 



we can now summarize the dominant contribution of the extended Higgs sector to the the 
coefficients a and b given by eqs. (H,^) [0: 



1 fB,/ 


' e Y 




167r2 2 \ 




1 fB, 




^ m 


167r2 2 




Ml. 



tl 



Y[Xt) — ^^^^ tan p- 



A 



ti 



Ml, 



tan^/? 



logr 



(24) 



where r = 1/yt = {Mh+ /rritY . Since > 1, the CP-odd neutral Higgs exchange interferes 
destructively with the SM contribution. Figures |l| show the contribution of the extended Higgs 
sector of the MSSM (assuming that sparticles contribute negligibly) or of the 2HDM(II)Q to 
BR{B^^ — > fJ^^fJ^^) as a function of Mh+ for different values of tan (3. These result which should 
be compared with the SM results 4 x 10~^ and 1 x 10~^°, respectively, agree for tan/5 > 30 with 



the ones given in and, for smaller values of tan /5, update the computations done earlier in 
refs. 



4. Chargino contribution 

Another source of amplification of the fiavour changing transitions induced by the CKM 
matrix is the chargino sector of the MSSM. Assuming that the squark mass matrices are diagonal 
in the super-CKM basis, the first result is that in the whole relevant parameter space the box 
diagram contribution to any of the Wilson coefficients remains small compared to the SM 
contribution. Furthermore, the penguin can change the predicted -Bi?(-B°^ — > /U^/i+) by no 
more than ~5-10% for tan/5 ~ 2 and ~20% for tan/5 ~ 0.5. The magnitude and sign of this 



^The well-known large radiative corrections to /i° and masses do not spoil these relations. Moreover, these 
corrections do not affect the neutral Higgs boson penguin contributions because they always modify significantly 
only the mass of that Higgs boson which almost does not couple to the down- type quarks and charged leptons. 

^In the case of the 2HDM(II) the subleading in tan /? contributions of the genuine vertex corrections in the CP- 
even Higgs boson penguin may be different than in the MSSM because the dimensionfuU couplings H^H~H^{h^) 
differ in both models Still, unless these coulings in the 2HDM(H) are very large (and numerically very 

different from their MSSM counterparts) so as to enhance the otherwise subleading contribution, figs ^ should 
be fairly representative also for the 2HDM(II) results. 
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Figure 1: Contribution of the Higgs sector of the MSSM or 2HDM(II) to BR{Bl^^^ ^+) as 
a function of the charged Higgs boson mass for tan/3 = 0.5 (sohd hues), 2 (dashed), 25 (dotted) 
and 50 (dot-dashed). 

contribution depends, apart from the masses of the sparticles involved, also on the chargino 
composition and on the mixing angle of the top squarks. For natural stop composition i.e. 
when the lighter stop is predominantly right-handed and the mixing angle is not too large 



1^], the chargino loop contribution to the Wilson coeffcients has opposite sign to that of 
the top quark loop and, hence, decreases the rate of the — decay. This is very 
similar to the opposite, as compared to the SM, sign of the chargino-stop loop contribution to 
Rh = r(Z° — > hh) /r{Z^ hadr) fl^ since, in view of the smallness of the box contribution, the 
two calulations are very similar We conclude that in the whole range of the MSSM parameter 
space the box and Z*^-penguin diagrams arising from chargino exchanges do not change the 
order of magnitude of the — > decay rate. 

Huge contribution to this rate can be however induced for large tan/3 > 30 by neutral 
Higgs boson penguins. This has been first maded explicit in ref. |^ in the approach based 
on the effective lagrangian method. Our numerical results are obtained by full computation 
of all relevant Feynman diagrams. Here we present only the derivation of the approximate 
formulae summarizing the dominant effects. To this end we consider the limit in which all 
soft SUSY breaking parameters, except for the ones which determine the Higgs potential, are 
much larger than the electroweak scale. In this limit, which allows us to work in the symmetric 
phase of the theory (i.e. with Vi = 0) in which sfermions still have definite chirality, we can 
construct the effective theory by integrating out sparticles (but not the Higgs fields). In this 
construction, threshold corrections shown in fig. ^ give rise to the effective Yukawa interactions 
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I 



d I 



'C 



*" — ' — z — — : — ' — * *" — ' — : — — : — ' — ■ 

ilu Hd dA H, d 

a) h) 

Figure 2: Diagrams giving rise to A^Yd i ^dYd respectively, in the construction of the effective 
theory. 



of the down-type quarks summarized by 

Ceff = - ei, {Yd + AdY.f^Hfqfd^'' - {A^Y.f^Hrqfd^'' + h.c. (25) 

where A, B are the generation indices and we work in the language of two-component Weyl 
spinors. In order to diagonalize the quark mass matrix arising after the electroweak symmetry 
breaking we perform first the standard CKM rotations (diagonalizing the original matrix Y^^) 
followed by the infinitesimal rotations 

{l + AVj^^f" ci^^^ci^«(l + AV^^y (26) 

with AVl \ AVff^ satisfying AVlr = -AVf^. Diagonal mass matrix for down-type quarks is 
obtained with 

- {A',Y,f^ + ^ {KY,f^ = {AV^Y^Yf + Y,\AVEr^ (27) 

where Y^ are already diagonal and A'^^j^-^Y^ are related to the original Au[d)Yd by the rotation 
diagonalizing the original Y^^. This leads [] to the effective Yukawa couplings of the neutral 
Higgs bosons of the form 

C = -^d' {-YdZ]^ + /^uYdZl' - tanPA'^YdZ},') d + h.c. 

+ ^d' {YdZ"^ + A'^YdZf + t^n^A'^^dZ]!) d H^^^ + h.c. (28) 

(In the lagrangian (|^) Y^ is diagonal and the rest of the notation is explained below eq. (|10|)) 
which in general generates the FCNC transitions. Note that the correction AdYd dissapeared 
altogether as it should, since it cannot cotribute to the FCNC transition. 

^Strictly speaking eq. ( p7|) must hold only for off diagonal elements; for A = _B the relation to^a = ^^d"""^ I 
is corrected but the net result is that in eq. ( p8| ) = ~\/2'm^^/vi again. 
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The correction A^y^^ in eq. (p5D is easily computed in the basis in which is diagonal 
and the initial Yj^^ = Y^V\q where Vab is the CKM matrix. Starting from the SUSY breaking 



part of the lagrangian |]T2 

. AB 



we obtain 



AB 



where Co is the standard three-point function 

1 



Co{a,b, c) 



a 



1 " 
log - - - 

a — c c 



b , b 
log- 

— c c 



(29) 



(30) 



Inserting (^) in eq. (|25|) and performing all steps making the usual assumption = 
Y^A^6^^ (i.e. that the trilinear soft terms are proportional to the Yukawa couplings) and 
keeping only the top Yukawa coupling leads to 



A' Y/^ 

u a 



1 



1671^ 



\tiY^Y^Atmc,C,{m'c^,Mi.Mi 



(31) 



where we have replaced /z with mc^ and the sign ± keeps track of the sign of /i. Using eq. (^ 
in eqs. (plIlOD yields the full vertex formfactors 



72k 



Ik , 



1 



-^KY,Z],Hanl3 



V2 



(32) 



(right formfactors are given by the Hermitean conjugation; they involve Yj and are, therefore, 
sub leading) . 

Detailed comparison of the above simplified calculation with the standard diagramatic ap- 
proach (in which one computes both, the self energy corrections the IPI vertex diagrams, in the 
phase in which the electroweak symmetry is broken) reveals that the dominant contributions 
given by eqs. (^) to the formfactors arise only from the self energy diagrams (the IPI vertex 
corrections contain one power of tan/3 less). Moreover, the comparison shows that one should 
replace At hj At = At + fj. cot P, Mi^ , Mi^ with the true mass eigenstates Mi^ , M^^ justifies 
the replacement of ±/i by the mass of the lighter chargino. 



Using eqs. 

fB 1 



■A 



) we get 
2m/ / e \ 



ti 



4 167r2 



Ml 



SwJ 

2mi 



Y{xt) - tan^ ± tan^ Mtm^.C 



W 



r - 1 Ml 



4 167r2'"*^M2, \sw) 



Ml o ^ log r Ml mf o ^ ~ 

^ tan^ /5— ^ ± TTT^TT^ tan^ (3Atmc,Co 



8M, 



w 



SM^Ml 



(33) 
(34) 
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Figure 3: BR{B^^^-^ — > /i /i"*") as a function of the stop mixing angle 6t for tan/? = 50, the 
hghter chargino mass 100 GeV and different values of Ma- Solid, dashed and dotted lines 
correspond to {Mf^, Mf^) equal (240,500), (400, 700) and (300, 850) GeV, respectively. In the 
left (right) pannels M2//X = 10(— 1), where M2 is the usual SU{2) gaugino mass parameter. 
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Figure 4: BR{B^f^^^ fJ'~fJ'^) as a function of the lighter chargino mass for tan/? = 50, 
Ma = 200 and 1000 GeV. Sohd (dashed) hues correspond to {Mf^,MiJ equal to {mc\,SmcJ 
and the stop mixing angle 9t — 10° (SO'') whereas dotted (dash-dotted) hnes to (3mci,5mcJ 
and 0t = 10° (30°), respectively. In the left (right) pannels = lO(-l). 
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Knowing that Y{xt) ~ 1 these formulae allow for a quick estimate of the effects. It is important 
to note that the contribution of chaginos to the rate grows as tan® (3 and therefore can be much 
larger than the contribution of the Higgs sector. 

Figure § shows the dependence of the full branching ratio BR{Bg —>■ including 
the SM, Higgs boson and chargino contributions, on the mixing angle of the top squarks for 
some values of the other MSSM parameters. The minimum around ~ corresponding to 
At ^ is clearly seen. Incidentally this plot also supports the replacement of /i by ±mc^ in 
eq. (pTD because very similarly (up to a reflection 6t — > —6t which follows from different signs 
of fi) looking curves in the left and right panels have the same but distinctly different /i. 
Another important feature of the chargino contribution is that it does not vanish if all sparticle 
mass parameters are scaled uniformly: M^. — AM^., mc. — XmCi, /i — > A/i, At XAt. This is 
clear from the fact that in such a case Cq X^^Cq- This is illustrated in figure ^ which shows 
BR{B^ ^^^) as a function of the lighter chargino mass for (M^^, MjJ equal to (mc, 3mc) 
and 6t = 10° (solid lines), (jnc,3mc) and 9t = 30° (dashed lines), (3mc,5mc) and 9t = 10° 
(dotted lines) and {3mc,5mc) and 9t = 30° (dash-dotted lines). In fact, keeping the stop 
mixing angle fixed requires that A^ scales as A^ rather than as A, which explains the growth of 
the rates with rrici ^g- S 

To check the correlation of the prediction for BR{B^^j^j f^ f^^) with the resits for BR{B —>■ 
Xs'~f) we have performed scans over the relevant parameter space of the MSSM. We took 
the following ranges: 100 < < 1000 GeV, 0.1 < iMz/yuj < 10, 1 < Mijmc, < 10, 

1 < M^JM^^ < 5 and —60° < 6t < 60° and rejected points for which ^psquarks > 6 x lO"'' 
and Mh < 107 GeV. For calculating BR{B — > Xg'j) we have used the routine based on refs. 



2T[| including the NLO matching conditions at the scale Mz for the top and charged Higgs 



contribution as in P^l ^^ and only the LO ones for the supersymmetric contribution 



We have not used the available NLO matching conditions for the supersymmetric particles since 
they are computed under the specific assumptions about the sparticle spectrum, not necessarily 
satisfied in the scan and, moreover, they are not valid for large values of tan f3. The theoretical 
uncertainty is taken into account by computing the rate for fih = 2.4 and 9.6 GeV and then by 
shifting its larger (smaller) value upward (downward) by the added in quadratures errors related 
to the uncertainties in as, mt, rric/mb, \VtbVt*/Vcb\^, and higher order electroweak corrections; 
we do not take into account the variation of the scale nw For a given set of the parameters of 
the MSSM the BR{B —>■ Xgj) value shown in fig. ^ corresponds to the lowest (highest) edge 
of the resulting band of theoretical predictions, if the whole band is above (below) the range 



allowed by CLEO |10|, and to the central point of the overlap of the theoretical and CLEO 



bands in the case such an overlap exists. 

The results of the scans, shown in fig. ^ demonstrate that the the CLEO result for BR{B 
Xs'j) does not ehminate the points corresponding to the largest values of BR^B'^f^^^ 1^1^^) 
and even does not exhibit any definite correlation between the two rates, especially for those 
points for which BR{Bg /W~/U~^) is very large. This is mainly due to the fact that the (LO) 
chargino contribution to BR{B —* Xs'j) decreases with growing sparticle masses whereas its 
contribution to BR{B^(^^^ f^ f^^) does not. f\ This allows to hope that even the full NLO 

®For the same reason we do not expect strong constraints from the K'^K^ or B^B^ mass differences, which 
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Figure 5: ^ X,7) versus BR{Bl^^^ for tan/? = 50, = 200 GeV in panels 

a) and b), tan/? = 50, Ma = 600 GeV in panel c) and tan/5 = 30, Ma = 200 in panel d). 
Limits from CLEO on BR{B X^'y) and on BR{B^ /^ A*^) are also shown by solid lines. 

computation of BR{B Xg'j) will not change this picture qualitatively. Moreover, fig. ^ shows 
that the present CLEO bound BR{B^ — >• ijT jX^) < 6.2 x 10~^ (shown in the upper-right plot 
by the vertical solid line) already puts some weak constraints on the MSSM parameter space 
in the case of large tan/? ~ rrit/mh and Ma ^ 300 GeV. 

5. Flavour changing induced by sfermion mass matrices 

Up to now we have assumed that the fermion and sfermion mass matrices are flavour 
diagonal in the same basis (the so-called super-CKM basis). In this section we consider the 
effects of nondiagonal entries in the sfermion mass matrices. It is customary to parametrize 
such nondiagonal entries by the so-called dimensionless mass insertions 0, 



- ■ (35) 

'■KJYY 



a priori also depend on the chargino and stop parameters. 
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where X,Y = L,R, K = u, d, I, {■M'k)xx diagonal elements of the XX blocks of the 

full mass squared matrices, and {AAij^)xY ^i-re the off diagonal entries of the XY blocks. Most 
of these insertions are bounded by the existing experimental data (for review, see refs. @, 0]). 

In the case of the l^l^ decay the relevant insertions are (^xy) ^"^^ (^xy) ' ^ ~ 

The first interesting point is to check the effects of the slepton mass insertions which are 
the only source of the decays — > //' (through the box diagrams with charginos in the loop) 

Very strong bounds from nonobservation of the transition /i — ^ 67 exist only on the ((^i^j) 
insertions whereas in the case of the B^ IV decay most important are the insertions 
5^j^j^ on which the bounds are weaker 0. Taking mc^ = 100 GeV, light stop Mf^ ^ 100 GeV 
and adjustnig the slepton sector mass parameters so to keep Mj > 90 GeV, M;> > 50 for 

/ 7 \ 13(23) 

[5[l) ~ 0.9 we get 

5i?(5° ^ [/') < 1.6 X 10'^^ 

BR{B°a II') < 3.8 X 10"^=^ (36) 

where //' = er or /xr. (The largest rates are obtained for |M2/yu| ^ 1 and small stop mixing 
angle 9t; the result scale approximately as l^i^^p.) For other parameters (heavier stops and 
charginos) branching fractions for these processes are, of course, smaller. 

We now discuss the effects of the flavour non-diagonal mass insertions in the down-type 
squark mass matrix and return to the // finel states. The approximate formulae accounting for 

the effects of the insertions (Sxy) easily derived in the so-called mass insertion method 

in which filavour off diagonal elements of the sfermion mass squared matrices are treated 
as additional interactions. Usually, the linear appromation in (oxy ) is sufficient to account 

/ , \/3 

for the results obtained with the full diagramatic calculation. In the case of a nonzero (^xy ) 
insertion the dominant contribution is expected to come from the diagrams involving gluinos, 
due to their strong coupling, = y/Anas, to quarks and squarks. Q (This expectation is 
confirmed by the numerical computation in which all one-loop contributions are taken into 
account.) Since at one-loop there are no box diagrams with gluinos we are left only with the 
and neutral Higgs boson penguin diagrams. As previously, the latter type of penguin diagrams 
is important only for large values of tan/5. Another important remark is that because the 
change of flavour in the gluino diagrams does not originate from the CKM mass matrix, the 
rates of the 5^ decays need not be suppressed compared to the rates of the B^ decays. 

For tan/5 values not too large, only the penguin contribution can be important. Di- 
rect computation shows however, that in the formula (0) terms linear in the mass insertions 

{^LL{RR)} cancel out completely between the self energy Ti^ and the proper vertex correc- 

9Ref. i gives (5^)'' < 0.2 {j^f, (^i^)'' < 700 and {S^^.f' < 100 {^f, which in most 

cases are superseded by (5^^ ^ 1 - the hmit in which one of the sneutrinos becomes tachionic. 

lOpor nonzero {Sxy)^^ insertion also neutrahnos contribute; moreover a nonzero {S'li^)^^ insertion induces, via 
the CKM matrix (see e.g. [|, H), nonzero {Sf^j^) insertions which affect in principle the chargino contribution. 
Both these effects are small but are taken into account in our numerical code. 
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tion Fx {X = L,R). Because of that, the effects of the nonzero {S'[^J^J^■^^ mass inser- 
tions, even taking into account their quadratic and higher contributions in gluino exchanges 
as well as neutralino diagrams, are small for tan/? values for which the neutral Higgs boson 

penguin graphs are negligible. Larger effects could come only from nonzero [^lrj mass 

insertions which, however, are strongly constrained P, | (^l_r) I < 0.07(mmaa;/lTeV), 

I (^Li?) I ^ 0.03(m^aa;/lTeV) (where rrimax = {max{Msq,mg)). Respecting these constraints, 
^ iBR{B°i fi-fi+) ) remains of order 4 x IQ-^ (10-^°). 

In the case of large tan (3 we have to compute in the linear approximation in the mass inser- 
tions both the scalar parts of the self energies and the IPI vertex corrections to the couplings 
djdjS'^^P'^). For the self energies we get 



- [AMiy^^mt {Ab + /itan/3) Do(m?, M^, M^, M|,)| (37) 

where Dq is the standard four-point function 

Do{a, b, c, d) = [Co{a, c, d) - Co{b, c, d)] , (38) 

a — 

rrig is the gluino mass and Md is the average mass of the two botom squarks. Similar formula 
is obtained for S;| with the replacement (AA^|))^^ — {AJ\4jy)j^j^. 

In the same approximation, for the vertex correction djdjP^ we get 

+ 7^T^(rn,,(AMiy'+(AMl)''mAzitanPDo{mlM (39) 



Zh' tan/?^-^A (m,, {^Ml)y^ + m,, (AAI?,) "J Doim% M^, M\ 



where for the three-linear soft term we have used = Y^A^. is similar, with (AA^|,)^'^ 
- {AMy)^R. Combining (|9|) and according to (|T3D, we see that (AA1|,)^^ cancels out. 
Moreover, since the CP-odd scalar A^ whose coupling to leptons is enhanced, corresponds to 
k = 1 and = cos/5 ~ 0, the second line in (|39|) is suppressed compared to the third one. 
Therefore, we can write: 



Fl ^ L\2^i!^tan2/?/i fAA^^)'"" mMml M^, M],, M^) (40) 



1 8 2 e rrih ^ 2 n f \ K j2\^-^ 

m^?,^'2sw Mw 

where we have retained only (AA^|,)^'^ which in F[ is multiplied by rridj = and neglected 
{AM."]^)^^^ which is multiplied by (in F^ it is the other way around). Similar calculation 
leads to 

" 1^^^'2^5^ ^'^D,{ml Ml, Ml, Ml) (41) 
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Figure 6: BR{B^^ for tan/? = 50, Ma = 200 GcV and (df^)^^ = 0.1 as a function of 

the n parameter. In the left panel {mg,At = Af,) equals: (300,0) GeV (solid line), (300,250) 
GeV (dashed), (800,0) GeV (dotted) and (800,250) GeV (dash-dotted); (m|)33 = (500 GeV)^, 
{ml)s3 = (m^)33 = (300 GeVf, {m\)KK = (600 GeV)^ for K ^ 3. In the right panel 
lmg,At = Ab) equals: (800,0) GeV (solid line), (800,450) GcV (dashed), (1500,0) GcV (dotted) 
and (1500,450) GcV (dash-dotted); (771^)33 = (900 GeV)^, (771^)33 = {mDsa = (700 GeV)^, 
{ml)KK = 1000 GeV2 for K 3. 



Computing the relevant Wilson coefficients we finally find for the coefficients a and b: 



a — 
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(42) 
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(43) 



in which we have displayed also the SM contribution to allow for easy estimate of the magnitude 
of the gluino contribution. It is essential that the dominant effect is due to the LL insertion and 
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not the LR one which is much more strongly constrained [Q, |^. Similarly as in the case of the 
chargino contribution through the neutral Higgs boson penguin graphs, also the gluino (and 
neutralino) contribution is proportional to tan® (3 and does not vanish when all SUSY mass 
parameters are uniformly scaled up (provided the dimensionless mass insertion is kept fixed). 
Figure |] shows the result of the full diagramatic computation of the SM and gluino exchange 
contributions to BR{Bg fJ'~fJ'~^) as a function of the /i parameter for (5^^)^^ = 0.1 and 
tan/3 = 50, Mj^ = 200 GeV. The minimum for /x = is clearly seen. The gluino contribution 
scales approximately as |((5ii)^^p. 

Figure ^ shows the results of the general scan over the MSSM parameter space in the 
form of the scatter plot BR{B Xgj) versus BR^B'^ /^~/^^) for (Sfj^)^^ = 0.1 and versus 
versus BR{B^ —* for (5^^)^^ = 0.05. The parameters have been varied in the following 

ranges: 100 < mc, < 600 GeV, 0.1 < |M2//i| < 5, m~g = SMa, -60° < Ot < 60°, Ab = At, 
0.5 < Mijmc, < 1.5, 1 < MijMi^ < 5, 0.25 < (m|))33/m? < 2.25. For other entries of 
the squark mass matrices we took {rn?x)KK = ("T'|))33- All points for which Mh < 107 GeV, 
^Psqark > 6 X lO"'^ (as wcll as poiuts with too light stops) have been rejected. Results for 
{^rr)^^ are similar. 



In agreement with the bounds given in refs. ||2|, ^, the measured by CLEO BR{B 
Xsj) does not constrain the rate of the BR{B^ A*"/^^) decay (nor does it exhibit any 
particular correlation with the latter) and the latter can attain values of the order of 10~^, 
respecting all the relevant phenomenological constraints. As expected, whenever the gluino 
contribution is dominant the rates of the B^ — > p~ and B^ — >■ fi~fi~^ decays are comparable 
which means that BR{B^ — > fj^~fi'^) can also be as large as 10~^ for {Sfj^Y^ = 0.1 (in the plot, 
we took (^ii)^^ = 0.05 in order to satisfy the bound (^^^(Hfl))^^ < 0.2(m^„^/lTeV) |, i for 
almost all points in the scan; however the biggest effects are for rrimax large for which (^^^(/jk))^'^ 
can be larger). It follows, that for such values of the MSSM parameters, the current CLEO 
bound, BR{B^ A^'/U^) < 6.2 x 10^^ |jTll puts constraints on {5'Il{rr)Y^ which are much 



stronger than the ones given in 
6. Conclusions 

We have performed a complete one loop diagramatic calculation of the decay rate of the -B^j-^., 
mesons into charged leptons. Both possible sources of the FCNC processe, the CKM mixing 
matrix and the off diagonal entries of the sfermion mass matrices, have been considered. For 
values of tan/?, for which the neutral Higgs boson penguin graphs are negligible the rates of 
these decays in the MSSM remain of the order of the SM prediction. 

Large enhancement of the SM prediction can occur for tan/3 ^ 1 provided the additional 
Higgs bosons predicted by the MSSM are not too heavy (all the large contributions behave as 
1/M\, where Ma is the mass of the CP-odd neutral Higgs boson). The contribution of the Higgs 
sector grows like tan^/3 and can give BR^B^^^ —>■ fi^fj^^) ~ 2x 10^^ for tan/3 ~ mt/mi,. Dominant 
effects of the chargino sector grow as tan® f3 and depend strongly on the top squark mixing. For 
tan/5 ~ mt/rrib and substantial mixing of the top squarks they can give BR{B^i^^-^ — >• 
up to 5 X 10~^(10~®) respecting other phenomenological constraints including the maesurement 
of BR[B Xs'j). Large effects, growing as tan®/? and exhibiting strong dependence on 
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Figure 7: BR{B X,7) versus ^ (with (5^^)^^ = 0-1) and BR{B^ fi-n+) 

(with (5ii)^^ = 0.05) for tan P = 50, Ma = 200 GeV in panels a) and b), tan P = 50, Ma = 600 
GeV in panel c) and tan /3 = 30, Ma = 200 in panel d). Limits from CLEO on BR{B X^-f) 
and on BR{B^ A*"/^"^) are also shown by solid hnes. 
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the /i parameter, can be also induced by the off diagonal elements of the down-type squark 
mass matrix. As we have shown, BR{B^(^^^ —>■ fi~fi~^) is sensitive to the 23 (13) off-diagonal 
entries of the LL and RR blocks of these matrices, which are not so strongly constrained 
by BR{B — > Xs7). For tan/? ~ mt/rrib and Ma ~ 200 GeV these effects can easily give 
BR{Bg /U^) larger than 10~^. It is also interesting that even for BR{B^ A* A*^) these 

effects can be so large that they could exceed the present CLEO limit which, therefore, 
already now puts constraints on the MSSM parameter space. 

Finally it is important to stress that, both types of effects growing as tan^/3 do not neces- 
sarily decrease as sparticles become heavy. However, they are sensitive to the mass scale of the 
extended Higgs sector. Thus, large deviation from the SM prediction observed in these decays, 
apart from being a signal of supersymmetry, would have important implications on the Higgs 
search at the LHC 
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